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Abstract 

The flux of low energy neutrinos {Ey ~ 4 GeV) has been studied with the MACRO 
detector at Gran Sasso via the detection of interactions inside the apparatus, and of 
upward-going stopping muons. Data collected in ~ 3 y with the full apparatus were ana¬ 
lyzed. The results are compatible with a deficit of the flux of atmospheric from below, 
and no reduction from above, with respect to Monte Carlo predictions. The deficit and 
the angular distributions are interpreted in terms of neutrino oscillations, and compared 
with the MACRO results on the upward throughgoing muons {Ey ~ 100 GeV) . 


Introduction 

Recent results from Super-Kamiokande (Fukuda 98) confirmed the anomaly in the 
ratio of contained to interactions, suggesting best £t parameters of sin^ 29 ~ 1.0 
and Am^ in the range of a few times 10“^ eV^ for disappearance. Other positive 
results come from Kamiokande (Fukuda 94), IMB (Casper 91) and Soudan 2 (Allison 
97) while earlier results from Frejus (Daum 95) and NUSEX (Aglietta 1989) are 
consistent with the expected number of contained events with smaller statistics. 
Also the results from MACRO (Ahlen 95) on up throughgoing muons presented an 
anomaly. 

The MACRO detector (Ahlen 93) is a large rectangular box (76.6 m x 12 m x 9.3 m) 
whose active detection elements are planes of streamer tubes for tracking and liquid 
scintillation counters for fast timing. The lower half of the detector is hlled with 
trays of crushed rock absorber alternating with streamer tube planes, while the 
upper part is open. The low energy flux can be studied by the detection of 

interactions inside the apparatus, and by the detection of upward going muons 
produced in the rock below the detector and stopping inside the detector (Fig. |^a). 
Because of the MACRO geometry, muons induced by neutrinos with the interaction 
vertex inside the apparatus can be tagged with time-of-flight {T.o.F.) measurement 
only for upgoing muons [IUpi= Internal Upgoing p). The downgoing muons with 
vertex in MACRO {IDp=Internal Downgoing p) and upward going muons stopping 
inside the lower part of the detector {UGSp = Upward Going Stopping p) can be 
identihed via topological constraints. Fig. lb shows the parent neutrino energy 
distribution for the three event topologies that can be detected by MACRO. 
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Internal Upgoing Events (lU) 

The data sample used for the Internal Upgoing (lU) events corresponds to a live-time 
of 3.16 years from April, 1994 up to November, 1997. During this period, ~ 22 10® 
atmospheric single muons were collected. The identihcation of lU/r events was based 
both on topological criteria and T.o.F. measurements. The basic requirement is the 
presence of at least two scintillator clusters in the upper part of the apparatus (see 
Fig. ^a) matching a streamer tube track reconstructed in space. A similar request 
was made in the analysis for the up throughgoing events (Ambrosio 98). 



Figure 1: (a) Sketch of different event topologies induced by neutrino interactions in 
or around MACRO. lU Internal Upgoing p; ID^= Internal Downgoing p; UGS^= 
Upgoing Stopping p; Up throughgoing = upward throughgoing p. In the figure, the stars 
represent the streamer tube hits, and the black boxes the scintillator hits. The time-of- 
flight of the particle can be measured for the lU^ and up throughgoing events, (b) Parent 
neutrino energy distributions for the three samples illustrated in (a). 

For lUp candidates, the track starting point must be inside the apparatus. To re¬ 
ject fake semi-contained events entering from a detector crack, the extrapolation of 
the muon track in the lower part of the detector must cross and not hre a mini¬ 
mum number of streamer tube planes and scintillator counters depending on track 
conhguration. The above conditions, tuned with Monte Carlo simulated events, 
account for detector inefficiencies and reduce the background contribution from up¬ 
ward throughgoing muons, which appear like semi-contained event to ~ 1%. The 
measured l/jS distribution is shown in Fig. ^ The measured muon velocity (3c is 
calculated with the convention that downgoing muons have 1/(3 values centered at 
-|-1 while upgoing muons have 1//3 at -1. It was evaluated that three events are 
due to an uncorrelated background. After background subtraction, 85 events are 
classihed as lU events. 
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Figure 2: The 1//3 distribution of the 85 lU events (dashed area). The ~ 1.2 10® events 
centered at 1//9 = +1 are downgoing atmospheric stopping muons. 

UpGoing Stopping (UGS) and Internal Downgoing (ID) Events 

The UGS+ID events were identified via topological constraints, and not with the 
T.o.F. For this analysis, the effective live-time is 2.81 y. The main request for the 
event selection is the presence of one reconstructed track crossing the bottom layer of 
the scintillation counters (see Fig. |l|a). All the hits along the track must be at least 
one meter inside from each wall of a MACRO supermodule. The selection conditions 
for the event vertex (or y stop point) to be inside in the detector are similar to those 
used for the lU search; they reduce to a negligible level the probability that an 
atmospheric muon produces a background event. The main difference with respect 
to the lU analysis (apart the T.o.F.) is that on average fewer streamer tube hits are 
present. To reject ambiguous and/or wrongly tracked events which passed the event 
selection, a systematic scan with the MACRO Event Display was performed. All 
real and simulated events which passed the event selection criteria were randomly 
merged. The finally accepted events had to pass the double scan procedure (the 
differences are included in the systematic uncertainty). In real events, three different 
subsamples were considered, according to the number of streamer tube hits, =NWP 
along the reconstructed track. 



E, 

(GeV) 

Rmin 
{g cm~‘^) 

% 

Events 

Bek. 

(±20%) 

MC 

(no osc.) 

A (All events) 

3.1 

30 

81 

241 

22 

256 

B {NWP > 3) 

3.9 

100 

90 

125 

5 

159 

C {NWP > 4) 

4.9 

160 

94 

66 

1 

95 


Table 1: Summary results (see text) for the three subsamples of ID-I-UGS events. 

As shown in Table 1, each subsample is distinguished by a different average parent 
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neutrino energy (col. 2 ); minimum range {g cm~‘^) of detector material to be crossed 
by the neutrino-induced particle (col. 3); percentage of C.C. interactions (col. 
4). In column 5 are given the numbers of real events. 

The main background source is due to upward going charged particles (mainly pi¬ 
ous) induced by interactions of atmospheric muons in the rock around the detector 
(Ambrosio 98). The background affects mainly sample A, and becomes smallest in 
sample C, as shown in col. 6 of Table 1. Sample B is used in the next section. 


Comparisons between Data and Monte Carlo 

The expected rates of lU and ID-I-UGS events were evaluated with a full Monte 
Carlo simulation. The and were allowed to interact in a large volume of 
rock containing the experimental Hall B and the detector. The rock and detector 
mass in the generation volume is 175 kton. The atmospheric flux given by the 
Bartol group (Agrawal 96) and the v cross sections of (Lipari 94) were used. The 
detector response was simulated using GEANT, and simulated events are processed 
in the same analysis chain as the real data. In the simulation, the parameters of the 
streamer tube and scintillator systems have been chosen in order to reproduce the 
real average efficiencies. The total theoretical uncertainty for the muon production 
from at these energies is ~ 25%. The systematic error is of the order of 10%, 
arising mainly from the simulation of detector response, data taking conditions, 
analysis algorithm efficiency, mass and acceptance of the detector. With our full 
MG, the prediction for lU events is 144 ± 36theor ± ^‘^syst, while the observed events 
are 85±9stat- The ratio R = {DATA/MC)iu = 0.59 ±0.06^tat iO.lbtw ± 0 . 06 ^ 3 ,^*. 
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Figure 3: Distribution of the cosine of the zenith angle (9) for ID-I-UGS (samples B) 
and lU events. The background-corrected data points (black points with error bars) are 
compared with the Monte Carlo expectation assuming no oscillation (full line) and two- 
flavour oscillation (dashed line) using maximum mixing and Am? = 2.5 10 “^ eU^. 


The MG prediction for ID-I-UGS events are shown in the last column of Table 
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1; in particular for sample B we expect 159 ± 4Qtheor ± ^^syst events. The ratio 
R = {DATA/MC)id+ugs = 0.75 ± 0.07stat ± 0.19^^ ± O-OSsyst- An almost equal 
number of ID and UGS neutrino induced events are expected in our data sample. 
The angular distribution of the lU and ID+UGS data samples, with the Monte 
Garlo predictions, are presented in Fig.|^. 

The low energy samples show a dehcit of the measured number of events in 
a uniform way over the whole angular distribution with respect to the predic¬ 
tions based on the absence of neutrino oscillations. The measured dehcit of low- 
energy events is in agreement with the MAGRO results on the throughgoing events 
(Ambrosio 98), i.e. with a model of disappearance with sin^20 ~ 1.0 and 
Am^ ~ 2.5 10“^ eV'^. In fact, the lU and UGS events have L ~ 13000 km, and in 
the energy range of few GeV the hux is reduced by a factor of two for maximum 
mixing and Am^ ~ 10“^ 10“^ eU^. No hux reduction is instead expected for ID 

events {L ~ 20 km). A global analysis of the diherent data sets is in progress. 
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